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ABSTRACT 

Q , We present cosmological hydrodynamical simulations of galaxy clusters aimed at studying 

V-j ' the process of metal enrichment of the intra-cluster medium (ICM). These simulations have 

C/3 I been performed by implementing a detailed model of chemical evolution in the Tree-SPH 

^J^. GADGET-2 code. This model allows us to follow the metal release from SNII, SNIa and AGB 

stars, by properly accounting for the lifetimes of stars of different mass, as well as to change 
the stellar initial mass function (IMF), the lifetime function and the stellar yields. As such, 
our implementation of chemical evolution represents a powerful instrument to follow the cos- 
mic history of metal production. The simulations presented here have been performed with 
CSJ . the twofold aim of checking numerical effects, as well as the impact of changing the model 

On ' of chemical evolution and the efficiency of stellar feedback. In general, we find that the dis- 

tribution of metals produced by SNII are more clumpy than for product of low-mass stars, 
If^ . as a consequence of the different time-scales over which they are released. Using a standard 

' Salpeter IMF produces a radial profile of Iron abundance which is in fairly good agreement 

, with observations available out to ~ 0.6i?5oo- This result holds almost independent of the 

■ numerical scheme adopted to distribute metals around star-forming regions. The mean age of 

enrichment of the ICM corresponds to redshift z ~ 0.5, which progressively increases outside 
the virial region. Increasing resolution, we improve the description of a diffuse high-redshift 

■ enrichment of the inter-galactic medium (IGM). This turns into a progressively more efficient 
H I enrichment of the cluster outskirts, while having a smaller impact at 0.5i?5oo- As for the 

effect of the model of chemical evolution, we find that changing the IMF has the strongest 
impact. Using an IMF, which is top-heavier than the Salpeter one, provides a larger Iron 
abundance, possibly in excess of the observed level, also significantly increasing the [0/Fe] 
relative abundance. Our simulations always show an excess of low-redshift star formation 
and, therefore, of the abundance of Oxygen in central cluster regions, at variance with ob- 
servations. This problem is not significantly ameliorated by increasing the efficiency of the 
stellar feedback. 

Key words: Cosmology: Theory - Galaxies: Intergalactic Medium - Methods: Numerical - 
X-Rays: Galaxies: Clusters 



1 INTRODUCTION 

Observations of the power spectrum of the CMB anisotropies, com- 
bined with those of the cosmic distance scales, and the statisti- 
cal properties of cosmic structures are now providing a convinc- 
ing validation of the standard cosmological scenario and placing 
quite stri ngent constraints on t he value of the cosmological param- 
eters (see lSpringel et al]|2006l . for a recent review, and references 
therein). As for the study of the formation process of cosmic struc- 



tures, this allows us to disentangle the effects of cosmic evolution 
from those induced by the astrophysical processes which determine 
the observational properties of galaxies and of the diffuse baryons. 
Star formation, feedback of energy and metals from explosions of 
supernovae (SN) and from Active Galactic Nuclei (AGN) play a 
fundamental role in determining the thermo-dynamical and chemo- 
dynamical status of cosmic baryons. 

In this respect, clusters of galaxies play a very important role. 
Thanks to the high density and temperature reached by the gas 
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trapped in their potential wells, they are the ideal signposts where to 
trace the past history of the inter-galactic medium (IGM). Observa- 
tions in the X-ray band with the Chandra and XMM-Newton satel- 
lites are providing invaluable information on the thermodynamical 
pro perties of t he intra-cluster medium (ICM; see lRosati et al .I2OO2I 
and lVoill2005l for reviews). These observations highlight that non- 
gravitational sources of energy, such as energy feedback from SN 
and AGN have played an important role in determining the ICM 
physical properties. 

At the same time, spatially resolved X-ray spectroscopy per- 
mits to measure the equivalent width of emission lines associated 
to transitions of heavily ionized elements and, therefore, to trace 
the pattern of chemical enrichment (e.g.. iMushotzkvl [20041 for a 
review). In turn, this information, as obtainable from X-ray ob- 
servations, are inextricably linked to t he history of formation and 
evolu tion of the galaxy population (e.g. lRenzinilll997| ; |PiDino et al.l 
I2OO2I , and references therein ), as inferred from obs ervational in 
the optical band. For instance, iDe Grandi et^ ( l2004h have shown 
that cool core clusters are characterized by a significant central en- 
hancement of the Iron abundance in the ICM, which closely corre- 
lates with the magnitude of the Brightest Cluster Galaxies (BCGs) 
and the temperature of the cluster. This demonstrates that a fully 
coherent description of the evolution of cosmic baryons in the con- 
densed stellar phase and in the diffuse hot phase requires properly 
accounting for the mechanisms of production and release of both 
energy and metals. 

In this framework, semi-analytical models of galaxy forma- 
tion provide a flexible tool to explore the space of parameters which 
describe a number of dynamical and astrophysical processes. In 
their most recent formulation, such models are coupled to dark 
matter (DM) cosmological simulations, to trace the merging his- 
tory of the halos where galaxy formation takes place, and include a 
treatment of metal producti on from type-la and type-II supemovae 
(SNIa and SNI I, hereafter; lOe Lucia et alj|2004 iNagashima et al.1 
l2005l :l Monaco et al. 200^, so as to properly add ress the study of 
the chemical enrichment of the ICM. ICoral 1 I2OO 6I recently applied 
an alternative approach, in which non-radiative SPH cluster simu- 
lations are used to trace at the same time the formation history of 
DM halos and the dynamics of the gas. In this approach, metals are 
produced by SAM galaxies and then suitably assigned to gas parti- 
cles, thereby pro v idin g a chemo-dynamical description of the ICM. 
iDomainko et al. I ( 12006) used hydrodynamical simulations, which 
include simple prescriptions for gas cooling, star formation and 
feedback, to address the specific role played by ram-pressure strip- 
ping in determining the distribution of metals. 

While these approaches offer obvious advantages with respect 
to standard semi-analytical models, they still do not provide a fully 
self-consistent picture, where chemical enrichment is the outcome 
of the process of star formation, associated to the cooling of the 
gas infalling in high density regions, as described in the numeri- 
cal hydrodynamical treatment. In this sense, a fully self-consistent 
approach requires that the simulations must include the processes 
of gas cooling, star formation and evolution, along with the corre- 
sponding feedback in energy and metals. 

A number of authors have presented hydrodynamical simula- 
tions for the formation of cosmic structures, which include treat- 
ments of the chemical evolution at different levels of complexity, 
using b oth Eulerian and S PH codes. Starting from the pioneering 
work bv lTheisetaLl ( ll992h . a number of chemo-dy namical models 
based on Eulerian codes have been presented (e.g., ISaniIancil998l : 
iRecchi et al]|200ll) with the aim of studying the metallicity evolu- 
tion of galaxies. Although not in a cosmological framework, these 



analyses generally include detailed models of chemical evolution, 

thereby accounting for the metal production from SNIa, SNII and 
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intermediate- and low-mass stars. Raiteri et al.l (1996) presented 

SPH simulations of the Galaxy, forming in a isolated halo, by fol- 
lowing Iron and Oxygen production from SNII a nd SNIa stars, 
also accounting for the effect of stellar lifetimes. iMosconi et al.l 
( 200 1 ) presented a detailed analysis of chemo-dynamical SPH sim- 
ulations, aimed at studying both numerical stability of the results 
and the enrichment prop erties of galactic objects in a cosmological 
context. iLia et alj ( I2OO2I) discussed a statistical approach to follow 
metal production in SPH simulations, which have a large number 
of star particles, showing appl ications to simulat i ons of a disc-like 
galaxy and of a galaxy cluster. lKawata & GibsonI ('2003) carried out 
cosmological chemo-dynamical simulations of elliptical galaxies, 
with an SPH code, by including the contri bution from SNIa and 
SNII, also accounting for ste llar lifet imes. Val daminil ( l2003h ap- 
plied the method presented bv lLia et"a l. (2002) to an exteded set of 
simulated galaxy clusters. This analysis showed th at profiles of the 
Iron a bundance are steeper than the observed ones. lTornatore et al] 
( I2OO4I) presented results from a first implem entation of a c hemi- 
cal evolution model in the GADGET-2 code ('Springel"2005'), also 
including the contribution from intermediate and low mass stars. 
Using an earlier version of the code presented in this paper, they 
studied the effect of changing the prescription for the stellar initial 
mass function (IMF) and of the feedback efficiency on the ICM 
enrichment in Iron, Oxygen and Silicon. A similar analysis has 
been presented by Romeo et al. (2006), who also considered the 
effect of varying the IMF an d the feedback efficiency on the en- 
richment pattern of the ICM. IScannapieco et al] ^2005 1) presented 
another implementation of a model of chemical enrichment in the 
GADGET-2 code, coupled to a self-consistent model for star for- 
mation and feedback (see also Scannapieco et al. ,2006) . In their 
model, which was applied to study the enrichment of galaxies, they 
included the contribution from SNIa and SNII, assuming that all 
long-lived stars die after a fixed delay time. 

In this paper we present in detail a novel implementa- 
tion of chemical evolution in the Tree-l-SPH GADGET-2 code 



dSpringel et id] I2OOII : ISpringell l2005b. which largely improves 



that originally introduced bv ^pringel & Hemguistl ( l2003ah (SH03 
hereafter). The model by SH03 assumes that feedback in energy 
and metals is provide d only by SNII , by assuming a Salpeter initial 
mass function (IMF; ISalpeteill955h , under the instantaneous recy- 
cling approximation (IRA; i.e. stars exploding at the same time of 
their formation). Furthermore, no detailed stellar yields are taken 
into account, so that the original code described a global metal- 
licity, without following the production of individual elements. Fi- 
nally, radiative losses of the gas are described by a cooling function 
computed at zero metallicity. 

As a first step to improve with respect to this description, we 
properly include life-times for stars of different masses, so as to 
fully account for the time-delay between star formation and release 
of energy and metals. Furthermore, we account for the contribution 
of SNII, SNIa and low and intermediate mass stars to the produc- 
tion of metals, while only SNII and SNIa contribute to energy feed- 
back. The contributions from different stars are consistently com- 
puted for any choice of the IMF. Also, radiative losses are com- 
puted by accounting for the dependence of the cooling function on 
the gas local metallicity. Accurate stellar yields are included so as 
to follow in detail the production of different heavy elements. The 
code implementation of chemical evolution is build in the GADGET- 
2 structure in an efficient way, so that its overhead in terms of com- 
putational cost is always very limited. 
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In the following of this paper, we will discuss in detail the 
effect that parameters related both to numerics and to the model 
of chemical evolution have on the pattern and history of the ICM 
chemical enrichment. While observational results on the ICM en- 
richment will be used to guide the eye, we will not perform here 
a detailed comparison with observations, based on a statistical en- 
semble of simulated clusters. In a forthcoming paper, we will com- 
pare different observational constraints on the ICM metal content 
with results from an extended set of simulated clusters. 

The plan of the paper is as follows. In Section 2 we describe 
the implementation of the models of chemical evolution. After pro- 
viding a description of the star formation algorithm and of its nu- 
merical implementation, we will discuss in detail the ingredients 
of the model of chemical evolution, finally reviewing the model 
of feedback through galactic winds, as implemented by SH03. In 
Section 3 we will present the results of our simulations. This pre- 
sentation will be divided in three parts. The first one will focus 
on the role played by numerical effects related to the prescription 
adopted to distribute metals around star forming regions (Sect. 3.1). 
The second part will concentrate on the numerical effects related to 
mass and force resolution (Sect. 3.2), while the third part (Sect. 3.3) 
describes in detail how changing IMF, yields, feedback strength 
and stellar life-times affects the resulting chemical enrichment of 
the ICM. The readers not interested in the discussion of the numer- 
ical issues can skip Sect. 3.1 and 3.2. We will critically discuss our 
results and draw our main conclusions in Section 4. 



2 THE SIMULATION CODE 

We use the Tr eePM-SPH GADGET-2 code jSpringel et"Zll200ll : 
ISpringell |2005|) as the starting point for our implementation of 
chemical evolution in cosmological hydrodynamical simulations. 
The GADGET-2 code contains a fully adaptive time-stepping, 
an explicit entropy-conserv ing formulation of the SPH scheme 
dSpringel & Hemauisl*2002l). heating from a uniform evolving UV 
background (Haardt & Madau 1996), radiative cooling from a zero 
metallicity cooling function, a sub-resolution model for star forma- 
tion in multi-phase interstellar medium (SH03), and a phenomeno- 
logical model for feedback from galactic ejecta powered by explo- 
sions of SNII. Chemical enrichment was originally implemented 
by accounting only for the con tribution of the SNII expected for 
a Salpeter IMF jSalpeteJI 19551 ). under the instantaneous recycling 
approximation (IRA) using global stellar yields. 

As we will describe in this Section, we have improved this 
simplified model along the following lines. 

(a) We include the contributions of SNIa, SNII and AGB stars to 
the chemical enrichment, while both SNIa and SNII contributes to 
thermal feedback. 

(b) We account for the age of different stellar populations, so that 
metals and energy are released over different time-scales by stars 
of different mass. 

(c) We allow for different Initial Mass Functions (IMFs), so as to 
check the effect of changing its shape both on the stellar popula- 
tions and on the properties of the diffuse gas. 

(d) Different choices for stellar yields from SNII, SNIa and PNe 
are considered. 

(e) Different schemes to distribute SN ejecta around star forming 
regions are considered, so as to check in detail the effect of chang- 
ing the numerical treatment of metal and energy spreading. 



2.1 The star formation model 

In the original GADGET-2 code, SH03 modeled the star formation 
process through an effective description of the inter-stellar medium 
(ISM). In this model, the ISM is described as an ambient hot gas 
containing cold clouds, which provide the reservoir of star forma- 
tion, the two phases being in pressure equilibrium. The density of 
the cold and of the hot phase represents an average over small re- 
gions of the ISM, within which individual molecular clouds cannot 
be resolved by simulations sampling cosmological volumes. 

In this description, baryons can exist in three phases: hot gas, 
clouds and stars. The mass fluxes between the these phases are reg- 
ulated by three physical processes: (1) hot gas cools and forms 
cold clouds through radiative cooling; (2) stars are formed from 
the clouds at a rate given a Schmidt Law; (3) stars explode, thereby 
restoring mass and energy to the hot phase, and evaporating clouds 
with an efficiency, which scales with the gas density. Under the as- 
sumption that the time-scale to reach equilibrium is much shorter 
than other timescales, the energy from SNe also sets the equilib- 
rium temperature of the hot gas in the star-formation regions. 

The original GADGET-2 code only accounts for energy from 
SNII, that are supposed to promptly explode, with no delay time 
from the star formation episode. Therefore, the specific energy 
available per unit mass of stars formed is esjve = gsnr x n'sN- 
Here, the energy produced by a single SN explosion is assumed to 
be esNe = 10^^ ergs, while the number of sta rs per solar mass 
ending in SNII for a Salpeter IMF ( ISalpeteJl95^, S55 hereafter) is 
u's'n = 0.0074 Mq\ 

In the effective model by SH03, a gas particle is flagged 
as star-forming whenever its density exceeds a given density- 
threshold, above which that particle is treated as multi-phase. 
Once the clouds evaporation efficiency and the star-formation 
(SF) timescale are specified, the value of the threshold is self- 
consistently computed by requiring (1) that the temperature of the 
hot phase at that density coincides with the temperature, Tu, at 
which thermal instability sets on, and (2) that the specific ejfec- 
tive energy (see eq. [11] of SH03) of the gas changes in a contin- 
uous way when crossing that threshold. Accordingly, the value of 
the density threshold for star formation depends on the value of the 
cooling function at Tu, on the characteristic time-scale for star for- 
mation, and on the energy budget from SNII. For reference, SH03 
computed this threshold to correspond to uh — 0.14 cm~"^ for the 
number density of hydrogen atoms in a gas of primordial composi- 
tion. 

In the simulations presented in this paper, we adopt the above 
effective model of star formation from a multi-phase medium. 
However, in implementing a more sophisticated model of chemi- 
cal evolution we want to account explicitly for stellar life-times, 
thereby avoiding the approximation of instantaneous recycling, as 
well as including the possibility to change the IMF and the yields. 
Therefore, while the general structure of the SH03 effective model 
is left unchanged, we have substantially modified it in several as- 
pects. Here below we describe the key features that we have imple- 
mented in the code, while we postpone further technical details to 
Sec. l [23] (. 

(1) The amount of metals and energy produced by each star par- 
ticle during the evolution are self-consistently computed for differ- 
ent choices of the IMF. In principle, the code also allows one to 
treat an IMF which changes with time and whose shape depends 
on the local conditions (e.g., metallicity) of the star-forming gas. 
This feature is not used in the simulations that we will discuss in 
this paper. 



4 Tornatore et al. 



(2) As in the original SH03 model, self-regulation of star forma- 
tion is achieved by assuming that the energy of short-living stars is 
promptly available, while all the other stars die according to their 
lifetimes. We define as short living all the stars with mass ^ Msl, 
where Ms l must be considered as a parameter whose value ranges 
from the minimum mass of core-collapse SNe (we assume 8M0), 
up to the maximum mass where the IMF is computed. This allows 
us to parametrize the feedback strength in the self-regulation of the 
star formation process. We emphasize that the above mass thresh- 
old for short living stars is only relevant for the energy available 
to self-regulate star formation, while metal production takes place 
by accounting for life-times, also for stars with mass ^ Msl - In 
the following runs we set Msl ~ SMq, thus assuming that all the 
energy from SNII is used for the self-regulation of star formation. 

(3) We include the contribution of metals to the cooling function. 
To this purpose, our implementation of cooling proceeds as follows. 
The original cooling function provided in the GADGET-2 code is 
used to account for the photo-ionization equilibrium of Hydrogen 
and Helium, while the tables by Sutherland & Dopita ( 1993) are 
used to account for the contribution of metals to th e cooling func- 
tion. We note that the cooling tables bv iSutherland & Dopit 3 j 19931) 
assume the relative proportions of different metal species to be so- 
lar. Including more refined cooling rates, which depend explicitly 
on the individual abundances of different metal species, involves a 
straightforward modification of the code. Due to the lack of an ex- 
plicit treatment of metal diffusion, a heavily enriched gas particle 
does not share its metal content with any neighbor metal poor par- 
ticle. This may cause a spurious noise in the cooling process, in the 
sense that close particles may have heavily different cooling rates, 
depending on how different is their metallicity. To overcome this 
problem, we decided to smooth gas metallicity using the same ker- 
nel used for the computation of the hydrodynamical forces (i.e., a 
B-spline kernel using 64 neighbors), but only for the purpose of the 
computation of the cooling function. Therefore, while each particle 
retains its metal mass, its cooling rate is computed by accounting 
also for the enrichment of the surrounding gas particles. 

(4) A self-consistent computation of the density threshold for star 
formation implies introducing a complex interplay between differ- 
ent ingredients. Firstly, changing the IMF changes the amount of 
energy available from short-living stars, in such a way that the 
threshold increases with this energy. Secondly, including the de- 
pendence of the cooling function on the local metallicity causes 
the density threshold to decrease for more enriched star form- 
ing gas. In the following, we fix the value of this threshold at 
hh ~ 0.1 cm^'^ in terms of the number density of hydrogen 
atoms, a value that has been adopted in a number of prev ious stud- 
ies of star formation in hydrodynamical simulations (e.g..lSummersl 
I1993I : lNavarro&Whitel ll993l ; iKatz et al.|[T993 : iKav et a l."2002'). 
and which is comparable to that, uh = 0.14 cm~^, computed by 
SH03 in their effective model. We defer to a future analysis the 
discussion of a fully self-consistent model to compute the star for- 
mation threshold. 

(5) The computation of the fraction of clouds proceeds exactly as 
in the original effective model (see eq.[18] in SH03), except that 
we explicitly include the dependence of the cooling function on the 
local gas metallicity, and the SN energy used for the computation 
of the pressure equilibrium is consistently computed for a generic 
IMF, as described above. 



2.2 The numerical implementation of star formation 

In order to the define the rule to transform star-forming gas par- 
ticles into collisionless star particles we closely follow the im- 
plementation by SH03 of the algorithm originally developed by 
Katz et al. ( 1996). This algorithm describes the formation of star 
particles as a stochastic process, rather than as a "smooth" contin- 
uous process. Basically, at a given time the star formation rate of 
a multi-phase gas particle is computed using a Schmidt-type law 
( ISchmidtll95"9h : 



(1) 



Here, x is the fraction of gas in cold clouds, so that xm is the mass 
of cold clouds providing the reservoir for star formation. Within the 
effective star formation model by SH03, the star formation time- 
scale, t*(p), is computed as 



(2) 



where is a parameter of the model, while pth is the density 
threshold for star formation, that we defined above. SH03 showed 
that the value of should be chosen so as to reproduce the 
observed relation between the disc-av eraged star form ation per 
uni t area and the gas surface de nsity jKennicut^ Il998l) . Follow- 
ing S pringel & Hernguistl ( l2003b!) . we assume tg = 1.5 Gyr, and 
checked that with this value we reproduce the Kennicutt law within 
the observational uncertainties. 

Following eq.l[Tll, the stellar mass expected to form in a given 
time interval At is 



= m 



(3) 



Within the stochastic approach to star formation, we define the 
number of stellar generations, TV,, as the number of star particles, 
which are generated by a single gas particles. Therefore, each star 
particle will be created with mass 



(4) 



where mg,o is the initial mass of the gas particles. Within this ap- 
proach, a star particle is created once a random number drawn in 
the interval [0, 1] falls below the probability 



m,,o 



exp 



(5) 



After the occurrence of a large enough number of star formation 
events, the stochastic star formation history will converge to the 
continuous one. In case a gas particle already spawned (A'^ — 1) 
star particles, then it is entirely converted into a star particle. As we 
shall discuss in Section Sec. l |2.3.U , star particles are allowed, in 
our implementation of chemical evolution, to restore part of their 
mass to the gas phase, as a consequence of stellar mass loss. Since 
this restored mass is assigned to surrounding gas particles, the latter 
have masses which can change during the evolution. Therefore, in 
the eq.lIU the actual mass rUg replaces the initial mass mgfi, which 
is assumed to be the same for all gas particles. As a consequence, 
star particles can have different masses due to (a) their mass loss, 
and (b) the mass of their parent gas particles. 

Clearly, the larger the number of generations, the closer is 
the stochastic representation to the continuous description of star 
formation. In the simulations presented in this paper we will use 
A*'* = 3 as a reference value for the number of stellar generations, 
while we checked that there is no appreciable variation of the final 
results when increasing it to TV, = 12. 
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2.3 The chemical evolution model 

Due to the stochastic representation of the star formation process, 
each star particle must be treated as a simple stellar population 
(SSP), i.e. as an ensemble of coeval stars having the same initial 
metallicity. Every star particle carries all the physical information 
(e.g. birth time t^, initial metallicity and mass), which are needed 
to calculate the evolution of the stellar populations, that they rep- 
resent, once the lifetime function (see Section l2.3.2l (. the IMF (see 
Section [23n ( and the yields (see Section [2331 ) for SNe and AGB 
stars have been specified. Therefore, we can compute for every star 
particle at any given time t > tb how many stars are dying as SNII 
and SNIa, and how many stars undergo the AGB phase, accord- 
ing to the equations of chemical evolution that we will discuss in 
Sec. ( 12.3. U below. The accuracy with which chemical evolution is 
followed is set by defining suitable "chemical" time-steps. These 
time-steps are adaptively computed during the evolution by fixing 
the percentage of SNe of each type, which explode within each time 
step. In our simulations, we choose this fraction to be 10 per cent 
for SNII and 2 per cent for SNIa. As such, these time-steps depend 
both on the choice of the IMF and on the life-time function. 

In the following, we assume that SNIa arises from stars be- 
longing to binary systems, having mass in the range 0.8-8 Mq 
dOreggio & Renzini 1983), while SNII arise from stars with mass 
> 8Mq. Besides SNe, which release energy and metals, we also 
account for the mass loss by the stars in the AGB phase. They con- 
tribute to metal production, but not to the energy feedback, and are 
identified with those stars, not turning into SNIa, in the mass range 
0.8-8 Mq. 

In summary, the main ingredients that define the model of 
chemical evolution, as we implemented in the code, are the follow- 
ing: (a) the SNe explosion rates, (b) the adopted lifetime function, 
(c) the adopted yields and (d) the IMF which fixes the number of 
stars of a given mass. We describe each of these ingredients here in 
the following. 

As we shall discuss in Sec. l |3.I| l, once produced by a star par- 
ticle, metals are spread over surrounding particles according to a 
suitable kernel. 



and Mbm, respectively. Fol lowing iGreggio & Renzinil h983h and 
iMatteucci & GreggioHl986l), we assume A = 0.1 , Msm = ZMq 
and Mbm = 16M(n.'Mat teucci & Gibsor] ( Il995l) applied a model 
of chemical enrichment of the ICM and found that A = Q.l was 
required to re produce the observed Iron enrichment, by assuming 
a Scalo IMF lscaldll986h . Changing the IMF would in principle 
require to change the value of A. Sin ce our model of ICM enrich- 
ment is quite different from that by Matteucci & GibsonI ( Il995h . 
we prefer here to fix the value of A and check the agreement with 
observations for different IMFs, rather than adjusting by hand its 
value case by case. Eq.Q holds under the assumption of impulsive 
star formation. Indeed, since each star particle is considered as a 
SSP, the associated star formation history, ip(t), is a delta-function, 
5(t — to), centered on the formation time to. 

As for the SNII and the low and intermediate mass stars stars, 
the rate is given by 



R 



SNII\AGB 



[t) = <t){m(t)) X 



dm{t) 
dt 



(V) 



We note that the above expression must be multiplied by a factor of 
(1 — A) for AGB rates if the interested mass m{t) falls in the same 
range of masses which is relevant for the secondary stars of SNIa 
binary systems. 

The release of energy and chemical elements by stars (binary 
systems in case of SNIa) of a given mass is obtained by multiplying 
the above rates by the yields pz. {m., Z), which give the mass of the 
element i produced by a star of mass m and initial metallicity Z. 
Then, the equation which describes the evolution of the mass pi (t) 
for the element i, holding for a generic form of the star formation 
history ip{t), reads: 



p,{t) = -i:{t)Z^{t) + 



^11 



(1 - A) ^^"'^ ^{t - r{m))pz, (m, ZMm) dm+ 



i)PZi (m, Z) dfi 



dm+ 



(8) 



2.3.1 The equations of chemical evolution 

We describe here the equations for the evolution of the rates of 
SNIa, SNII and AGB stars, along with their respective metal pro- 
duction. We provide here a short description of the basic results 
and of the equations, which are ac tually solved in o ur simulations, 
while we refer to the textbook by iMatteuccH ( |2003|) for a detailed 
discussion. 

Let r(m) be defined as the life-time function, i.e. the age at 
which a star of mass m dies. Accordingly, the the rate of explosions 
of SNIa reads 



RsNIa{t) = — 



dm(t) 



(6) 



m2— T l(t) 

Mbm 



24 mi A 



{ruB)- 



-dmi 



where the term in the first line is the mass of stars dying at the time 
t, T~^{t) is the inverse of the lifetime function r(m), (^(m) is the 
IMF and A is the fraction of stars in binary systems of that partic- 
ular type to be progenitors of SNIa. The integral is over the mass 
ma of the binary system, which runs in the range of the minimum 
and maximum allowed for the progenitor binary system, AlBm 



jMb^ _ T{m))pz^(m, Z)ip{m) dm+ 

^Mbm ^ '^i^))Pz, ("^1 Z)ip{m) dm. 

In the above equation. Ml and Mu are the minimum and maxi- 
mum mass of a star, respectively. In the following we use M l ~ 
O.IMq and Mu = WOMq. The term in the first fine of eq.lUll ac- 
counts for the metals which are locked up in stars. The term in the 
second line accounts for metal ejection contributed by SNIa. Here 
we have explicitly written the inner integral that accounts for all the 
possible mass ratios /j, — m2/{m\ -\- 7712) between the secondary 
star mass and the total mass; /imin is the minimum value of p and 
f{p) is the corresponding distribution function. The terms on the 
third and fourth lines describe the enrichment by mass-loss from 
intermediate and low mass stars, while the last line accounts for 
ejecta by SNII. 

The p distribution function is assumed to be 



= 2^+^(1 + 7)m 



(9) 



where 7 = 2. This functional form of f[p) has been derived from 
statistical studies of th e stellar po pulation in the solar ne ighborhood 
jTutukov & lungelson 1980; Ma tteucci & Recchil200ll) . The value 
of /imiii is calculated for a binary system of mass Mb as 
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fJ-n 



/ 7712 niB — 0.5Mbm 



(10) 



Taking the impulsive star-formation, the terms in eq.lts} must be 
recast in the form that we actually use for calculating the rates. 

In order to solve eqs.Q, JT]) and ^ in the GADGET-2 code we 
proceed as follows. At the beginning of each run, two tables, one 
for SNII and one for SNIa and low and intermediate mass stars, are 
built to specify at what delay times the chemical evolution should 
be calculated. The accuracy of these "chemical time-steps" is set 
by two run-time parameters that specify what fraction of stars must 
be evolved at each step. Accordingly, during each global time-step 
of the code only a small fraction (typically few percent) of all stars 
is processed. This implementation of chemical evolution is quite 
efficient in terms of computational cost, especially when the num- 
ber of stars grows. We verified that using A'^* = 3 for the number of 
stellar generations, the overhead associated to the chemical evolu- 
tion part amounts only to ^ 10 per cent of the total computational 
cost for a typical simulation. 



2.3.2 The lifetime function 

In our reference run, w e use the function given by 
IPadovani & Matteuccil h993h (PM hereafter). 



r(m) 



■|^q[(1. 34-^1. 79-0. 22(7. 76-log(m)))/0. 11] -9 

for m ^ 6.6 Mq 



(11) 



1.2 m" 



-I- 0.003 otherwise 



Furtherm ore, we also consid e r the l ifetime function originally pro- 
posed by iMaeder & Meynet' (198?) (MM hereafter), and extrapo- 
lated by [chiappini et al. ( 1997.) to very high (> 6OM0) and very 
low (< 1.3Af0) masses: 



Tim) = < 



10 
10 



10^ 



-3.71ogm+1.351 



-2.51 log m+0.77 



-1.78 log m+0.17 



m 1.3 M0 
1.3 < m ^ 3M0 



3 < m ^ 7 M, 







(12) 



7 < m < 15 Mg 



10-" **" m-0.94 15 < m < 53 M0 

1.2 X m-'^-^^ + 0.003 otherwise 

We refer to the paper bv lRomano et alj l l2005h for a detailed discus- 
sion on the effects of different lifetime functions on the chemical 
enrichment model of the Milky Way. 

A comparison between the life-time functions of egs. dl It and 
( I12t is shown in Figure [T] The main difference between these two 
functions concerns the life-time of low mass stars (< 8 Mq). The 
MM function delays the explosion of stars with mass ^ 1 Mq , 
while it anticipates the explosion of stars below 1 Mq with respect 
to PM function. Only for masses below IM0 , the PM function pre- 
dict much more long-living stars. We have verified that, assuming a 
Salpeter IMF (see below), the SNIa rate from a coeval stellar popu- 
lation is expected to be higher after ~ 3 Gyr when the MM lifetime 
function is adopted. This implies that different life-times will pro- 
duce different evolution of both absolute and relative abundances. 
This will be discussed in more detail in Sect. 13. 3.3] 

We point out that the above lifetime functions are indepen- 
dent of metallicity, whereas in principle this dependen ce can be in- 
cluded in a model of chemical evolution. For instance, iRaiteri et al.l 
d 19961) used the metallicity-d ependent lifetimes as obtained from 
the Padova evolutionary tracks iBertelli et ah 
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Figure 1. The dependence of the lifetime functions on the stellar mass. The 
solid and the dashed lines are for the lifetimes by IPadovani & Matteuccil 
a993i) and bv lMaeder & Mevnell jl989l) . respectively. 



2.3.3 Stellar yields 

The stellar yields specify the quantity pzi{m, Z), which appears 
in eq[8]and, therefore, the amount of different metal species which 
are released during the evolution of each star particle. In the 
runs presented in this work, we ad opt the yields provided by 
Ivan den Ho ek & Groenewegen ('1997') the low and intermediate 
mass stars and by Thielemann et al. (2003) fo r SNIa. As for SNII 
we ad opt the metallicity-dependent yields by IWooslev & Weaver! 
( Il995h in the reference ru n, while we will also show the effect of 
using instead the yields bv lChieffi & Limongil ( l2004h (WW and CL 
respectively, hereafter; see Table |2j. We also assume that all the 
stars having masses > 40 A'Iq directly end in black holes. 

Along with freshly produced elements, stars also release non- 
processed elements. Sometimes, papers presenting yields table do 
not explicitly take into account thes e non-processed yields (e.g., 
Ivan den Ho ek & Gro enewegenI 19971) . In order to account for them, 
whenever necessary, we assume that the non-processed metals are 
ejected along with Helium and non-processed Hydrogen. 

Besides H and He, in the simulations presented in this paper 
we trace the production of Fe, O, C, Si, Mg, S. The code can be 
easily modified to include other metal species. 

In Figure|2]we show the ratios between the abundances of dif- 
ferent elements, as expected for the WW and the CL yields, from 
the SNII of a SSP Different curves and symbols here correspond 
to different values of the initial metallicity of the SSP. Quite appar- 
ently, the two sets of yields provide significantly different metal 
masses, by an amount which can sensitively change with initial 
metallicity. In Sect |3.3.2| we will discuss the effect of changing 
yields on the resulting enrichment of the ICM and of the stellar 
population in simulated clusters. 



2.3.4 The initial mass function 

The initial mass function (IMF) is one of the most important quan- 
tity in modeling the star formation process. It directly determines 
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Figure 2. The ratio A/J^'^/Mj^^ between the mass of species j, pro- 
duced by the SNI I of a SSP, when using the two sets of yields by 
IWoosIev & WeaveJ j 19951) and by Chieffi & Limongi (2004) for different 
values of the initial SSP metallicity. Different symbols are for different val- 
ues of the initial metallicity of the SSP, as reported by the labels. 
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Figure 3. The dependence of the shape of the IMF on the stellar mass, rel- 
ative to the Salpeter IMF iSalpeter 1955). Dashed and dot-dashed curves 
are for the IMF by Arimoto & Yoshii ( 1987) and by Kroupa 1 2001 ), respec- 
tively. The horizontal solid line indicates for reference the Salpeter IMF. 



the relative ratio between SNII and SNIa and, therefore, the rela- 
tive abundance of a-elements and Fe-peak elements. The shape 
of the IMF also determines how many long-living stars will form 
with respect to massive short-living stars. In turn, this ratio affects 
the amount of energy released by SNe and the present luminosity 
of galaxies, which is dominated by low mass stars, and the (metal) 
mass-locking in the stellar phase. 

As of today, no general consensus has been reached on 
whether the IMF at a given time is universal or strongly dependent 
on the environment, or wheter it is time-dependent, i.e. whether lo- 
cal variations of the values of temperature, pressure and metallicity 
in star-forming regions affect the mass distribution of stars. 

Nowadays, there are growing evidences that the IMF in the 
local universe, expressed in number of stars per unit logarithmic 
mass interval, is likely to be a power-law for m* > 1 Mq with 
slope X ~ 1.35, while it becomes flat below the IMq thresh- 
old, possibly even taking a negative slope b elow ^0.2 (e.g., 
lKroupall200lh . Theoretical arguments (e.g., lLarson|[l998[) suggest 
that the present-day characteristic mass scale ~ 1 Mq should have 
been larger in the past, so that the IMF a t higher redshift was top- 
heavier than at present. IChiappini et"ai] [200Q) showed that vary- 
ing the IMF by decreasing the characteristic mass with time, leads 
to results at variance with observations of chemical properties of 
the Galaxy. While the shape of the IMF is determined by the lo- 
cal conditions of the inter-stellar medium, direct hydrodynamical 
simulations of star formation in molecular clouds are only now ap- 
proaching the required resolution an d sophistication level t o make 
credible predictions on the IMF (e.g., iBate & Bonneiill2005h . 

In order to explore how the IMF changes the pattern of metal 
enrichment, we implement it in the code in a very general way, so 
that we can easily use both single-slope and multi-slope IMFs, as 
well as time-evolving IMFs. In this work we use single-slope IMFs 
defined as 



using X — 1.35 for the standard Salpeter IMF ( ISalDetedll955h in 
our reference run. In the above equation, N is the number of stars 
per unit logarithmic mass interval. We will explore the effect of 
changing the IMF by als o using a top-heavier IMF with x = 0.95 
l l Arimoto & Yoshiill 19871 AY hereafter), as well as the multi-slope 
IMFbv lKroupal ( l200lh . which is defined as 



(m) oc 



-1.7 
-1.2 
-0.3 



m> IM 







m 

m '^'^ O.5s;m<lM0 



(14) 



m < 0.5 Mq 



(;/>(m) = d 7V/dlog m oc m 



(13) 



We show in Figure [3] the number of stars, as a function of 
their mass, predicted by the AY and Kroupa IMFs, relative to those 
predicted by the Salpeter IMF. As expected, the AY IMF predicts 
a larger number of high-mass stars and, correspondingly, a smaller 
number of low-mass stars, the crossover taking place at ~ 2Mq . 
As a result, we expect that the enrichment pattern of the AY IMF 
will be characterized by a higher abundance of those elements, like 
Oxygen, which are mostly produced by SNII. On the other hand, 
the Kroupa IMF is characterized by a relative deficit of high-mass 
stars and, correspondingly, a relatively lower enrichment in Oxygen 
is expected. 

Since clusters of galaxies basically behave like "closed 
boxes", the overall level of enrichment and the relative abundances 
should directly reflect the condition of star formation. While a num- 
ber of studies have been performed so far to infer the IMF shape 
from the enrichment patter n of the I CM, no general consensus has 
been reached. F or instance. lRenzinil ( ll997l) . |Pipino et al] ( l2002h and 
iRenzinil ( l2004b argued that both the global level of ICM enrich- 
ment and the a/Fe relative abundance can be reproduced by assum- 
ing a Salpeter IMF, as long as this relative abundance is subsolar 
Indee d, a different conclusion has been reached by other authors 
(e.g., iLoewenstein & Mushotzkvlll996h in t he attempt of explain- 
ing [a/Fe] > in the ICM. For instance. [Portinari et al] ( l2004h 
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used a phenomenological model to argue that a standard Salpeter 
IMF can not account for the observed g/Fe ratio in the I CM. A 
similar conclusion was also reached by Nagashim a et al ] i l2005h . 
who used semi-analytical models of ga laxy formation to tra ce the 
production of heavy elements, and by iRomeo et al. I ( l2006h , who 
used hydrodynam ical simulations including chemical enrichment. 
ISaroetal.l ( [2006h analysed the galaxy population from simulations 
similar to those presented here. This analysis led us to conclude that 
runs with a Salpeter IMF produce a color-magnitude relation that, 
with the exception of the BCGs, is in reasonable agreement with 
observations. On the contrary, the stronger enrichment provided by 
a top-heavier IMF turns into too red galaxy colors. 

In summary, our impleme ntation of chemic a l evol ution is 
quite simil a r to th at adopted bv lKawata&GibsorJ j2003h and by 
iKobavashil ( |2004|) . while diff erences exist with r espe ct to other 
implem entations . For instance. iRaiteri et al.l il99^ and lValdaminil 
l l2003 h also used a scheme similar to ours, but n eglected the con- 
tribution from low- and intermediate-mass stars. iLia et alj Bmh 
adopted a coarse stochastic description of the ejecta from star parti- 
cles: differently from our continuous description, in which enriched 
gas is continuously released by star particles, they assumed that 
each star particle is assigned a given probabilit y to be entire l y con- 
vert ed into an enriched gas pa rticle. Finally. Mosconi et alj i200lh 
and IScannapieco et al] ( l2005h neglected delay times for SNII, as- 
sumed a fixed delay time for SNIa and neglected the contribution 
to enrichment from low- and intermediate-mass stars. 



2.4 Feedback through galactic winds 

SH03 discussed the need to introduce an efficient mechanism to 
thermalize the SNe energy feedback, in order to regulate star for- 
mation, thereby preventing overcooling. To this purpose, they in- 
troduced a phenomenological description for galactic winds, which 
are triggered by the SNII energy release. We provide here a basic 
description of this scheme for galactic winds, while we refer to the 
SH03 paper for a more detailed description. The momentum and 
the kinetic energy carried by these winds are regulated by two pa- 
rameters. The first one specifies the wind mass loading according 
to the relation, Mw = f?A/t, where M, is the star formation rate. 
Following SH03, we assume in the following r] = 2>. The second 
parameter determines the fraction of SNe energy that powers the 
winds, ^MwVw = X^SNeMt,, where esNe is the energy feed- 
back provided by the SNe under the IRA for each A-Iq of stars 
formed. In the framework of the SH03 effective model for star for- 
mation, winds are uploaded with gas particles which are stochas- 
tically selected among the multiphase particles, with a probability 
proportional to their local star formation rate. As such, these par- 
ticles come from star-forming regions and, therefore, are heavely 
metal enriched. Furthermore, SH03 treated SPH particles that be- 
come part of the wind as temporarily decoupled from hydrodynam- 
ical interactions, in order to allow the wind particles to leave the 
dense interstellar medium without disrupting it. This decoupling is 
regulated by two parameters. The first parameter, pdcc, defines the 
minimum density the wind particles can reach before being coupled 
again. Following SH03, we assumed this density to be 0.5 times the 
threshold gas density for the onset of star formation. The second pa- 
rameter, Idcc, provides the maximum length that a wind particle can 
travel freely before becoming hydrodynamically coupled again. If 
this time has elapsed, the particle is coupled again, even if it has 
not yet reached pdec. We assumed Idcc = 10 h~^kpc. 

While we leave the scheme of kinetic feedback associated to 



wind galactic unchanged, we decide to use instead the value of the 
wind velocity, Vu,, as a parameter to be fixed. For the reference run, 
we assume ~ 500 km s^^. With the above choice for the wind 
mass loading and assuming that each SN provides an energy of 
10^^ ergs, this value of corresponds to esNe — 1 for a Salpeter 
IMF. We will also explore which is the effect of assuming instead 
a stronger feedback, with Vuj = 1000 km s~^, on the pattern of 
chemical enrichment. 



3 RESULTS 

In this Section we will discuss the results of a set of simula- 
tions of one single cluster. The cluster that we have simulated has 
been extracted from a low-resolution cosmological box for a flat 
ACDM cosmology with Qrn ~ 0.3 for the matter density parame- 
ter, h = 0.7 for the Hubble parameter, = 0.039 for the baryon 
density parameter and as — 0.8 for the normalization of the power 
spectrum of density pe rturbations. Using th e Zoomed Initial Con- 
dition (ZIC) technique jTormen et alJI 19971) , mass resolution is in- 
creased in the Lagrangian region surrounding the cluster, while also 
high-frequency modes of the primordial perturbation spectrum are 
added. The main characteristics of this cluster (Cll hereafter) are 
summarized in Table [T] along with the mass resolution and the 
Plummer-equivalent softening parameter used in the runs (note that 
the softenings are set fixed in physical units from z = 2 to 2 = 0, 
while they are fixed in comoving units at higher redshift). 

The reference run of Cll is performed under the follow- 
ing assumptions: (a) metals produced by a star particle are dis- 
tributed to surrounding gas particles by using a SPH spline kernel 
with density weighting over 64 neighbors for the distribution of 
metals around star forming p articles (see Sect. 13. U: (b) Salpeter 
IMF; (c) stellar yields fro m Thiele mann et al.l ( 12003 1 for SNIa, 
Ivan den Hoek & Groenewegen 1 1997i) for the low and intermediate 
mass stars an dlWooslev & WeaveJ (Il995h for SNII; (dj Life-time 
function from lPadovani & Matteuccif l 1993h ; (e) v,„ — 500 kms^^ 
for the velocity of winds. In the following we will show the ef- 
fect of changing each one of these assumptions. Therefore, we will 
explore neither the cluster-by-cluster variation in the enrichment 
pattern of the ICM, nor any dependence on the cluster mass. We 
perform instead a detailed study of the effect of changing a number 
of parameters which specify both the numerical details of the im- 
plementation of the model of chemical evolution model. We defer 
to a forthcoming paper the properties of the ICM enrichment for a 
statistical ensemble of simulated galaxy clusters. We summarize in 
Table[2]the aspects in which the various runs of Cll differ from the 
reference one. 

In addition, we also performed simulations of another cluster 
(C12 in Table[Tll, whose initial conditions are generated at three dif- 
ferent resolutions, for the same cosmologi cal model, by spa nning 
a factor 15 in mass resolution (see also B organi et al. 1 120061 . The 
lowest resolution at which the C12 cluster is simulated is compa- 
rable to that of the Cll runs. As such, this series of runs, that is 
performed for the same setting of the Cll reference (R) run, allows 
us to check in detail the effect of resolution on the ICM chemical 
enrichment. 

Besides global properties, we will describe the details of the 
ICM enrichment by showing radial profiles of the Iron abundance 
and of the relative abundances of [0/Fe] and [Si/Fefl the his- 



^ We follow the standard bracket notation for the relative abundance of 
elements A and B: [A/B] = \o^{ZaIZb) - \og{Z a,q / Zb,q). 
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Cluster M200 R200 



Cll 2.2 0.98 5.7 5.0 

C12 1.4 0.85 

Rl 2.31 5.2 

R2 0.69 3.5 

R3 0.15 2.1 



Table 1. Basic properties of the two simulated clusters Column 1 : cluster 
name; Column 2: mass within the radius, R200, encompassing an overden- 
sity of 200 times the critical density, pc ; Column 3 : value of R200 ; Column 
4: initial mass of the gas particles; Column 5: Plummer-equivalent soften- 
ing of the gravitational force aX z = 0. 



R 


Reference run 


N16/N128 


B-spline kernel for metal distribution with 16/128 neighbors 


MW 


B-spline kernel for metal distribution using mass weighting 


TH 


Top-hat kernel for metal distribution 


AY 


Top-heavier IMF (Arimoto & Yoshii 1987) 


Kr 


Kroupa IMF (Kroupa 2001) 


CL 


Yields for SNII from Chieffi & Limongi (2004) 


NW 


No feedback associated to winds 


sw 


Strong winds with = 1000 kms^^ 


MM 


Life-time function by Maeder & Meynet (1989) 



Table 2. Characteristics of the different CI 1 runs with respect to the refer- 
ence run (see text). Column 1; Name of the run; Column 2: characteristic of 
the run. 



tograms of the fraction of gas and stars having a given level of 
enrichment, and the time evolution of the metal production and en- 
richment. A comparison with observational data will only be per- 
formed for the abundance profiles of Iron, which is the element 
whose distribution in the ICM is best mapped. As a term of com- 
parison, we will use a subset o f 5 clusters, out of the 16, observed 
with Chandra and analysed bv lVikhUnin et alj j2005b . These clus- 
ters are those having temperature in the range 2-A keV, comparable 
to that of the simulated Cll and C12 clusters. Since the profiles of 
Fe abundance are compared with observations, we plot them out to 
^socd which is the maximum scale sampled by observations. We 
plot instead profiles of [Si/Fe] and [0/Fe] out to 2i?virin order to 
show the pattern of relative enrichment over a large enough range 
of scales, where the different star formation histories generates dif- 
ferent contributions from stars of different mass. Here and in the 
following we assume that a SPH particle belongs to the hot diffuse 
gas of the ICM when it meets the following three conditions: (i) 
whenever the particle is tagged as multiphase, its hot component 
should include at least 90 per cent of its mass; (ii) temperature and 
density must not be at the same time below 3 x 10*K and above 
500pbar respectively, where ptar is the mean cosmic baryon den- 
sity. While observational data on the Iron abundance profiles are 
used to "guide the eye", we emphasize that the primary aim of this 
paper is not that of a detailed comparison with observations. This 
would require a statistically significant set of simulated clusters, 
sampling a wider range of temperatures, and is deferred to a forth- 
coming paper. 



3.1 The effect of metal spreading 

A crucial issue in the numerical modeling of ICM enrichment con- 
cerns how metals are distributed around the star particles, where 
they are produced. In principle, the diffusion of metals around 
star forming regions should be fully determined by local physical 
processes, s uch as thermal motio ns of ions (e.g., Sarazin 1988), 
turbulence 1 Rebusco et al.l |2005|) , ram-p ressure stripping (e.g., 
Domainko et al l2006l ). galactic ejecta (e.g.. lStrickland & Heckmarj 
2006h . However, for these mechanisms to be correctly operating in 
simulations, it requires having high enough resolution for them to 
be reliably described. A typical example is represented by turbulent 



^ In the following we will indicate with the radius encompassing an 
overdensity of Apcr, where pcr is the critical cosmic density. In this way, 
will be defined as the mass contained within Ha- Furthermore, we 
will define /Jvirto be the radius encompassing the virial overdensity Avir, 
as predicted by the spherical top-hat collapse model. For the cosmology 
assumed in our simulations, it is Avir — 100. 



gas motions (e.g.. lRebusco et al]|2005l) . The description of turbu- 
lence requires not only sampling a wide enough dynamical range, 
where the energy cascade takes place, but also a good knowledge of 
the plasma viscosity, whi ch determine the sc ale at which turbulent 
energy is dissipated (e.g., iDolag et alj|2005h . While approaches to 
include physical viscosity in the SPH scheme have been recently 
pursued (Sijacki & Springel 2006b), a full accounting for its ef- 
fect requires including a self-consistent description of the magnetic 
field which suppresses the mean free path of ions. 

Owing to the present limitations in a fiilly self-consistent nu- 
merical description of such processes, we prefer to adopt here a 
simplified numerical scheme to describe the diffusion of metals 
away from star forming regions. We will then discuss the effect of 
modifying this scheme, so as to judge the robustness of the results 
against the uncertainties in the modeling of the metal transport and 
diffusion. 

Every time a star particle evolves, the stars of its population, 
the ejected heavy elements and energy must be distributed over the 
surrounding gas. The code accomplishes this task by searching for 
a number Ng of gas neighbours and then using a kernel to spread 
metals and energy according to the relative weights that the kernel 
evaluation assigns to each gas particle. In this way, the fraction of 
metals assigned to the i-th neighbour gas particle can be written as 



(15) 



where Wt is kernel value at the position of the i-th particle A^„i, 
is the number of neighbors over which metals are distributed and 
the denominator enforce conservation of the mass in metals. As for 
the weighting function W, the most natural choice, which is the 
one usually adopted in the SPH chemo-dynamical models so far 
presented in the literature, is to use the same B-spline kernel used 
for the computation of the hydrody namical quantities, a lso using 
the same number of neighbors (e.g., iMosconi et al.ll200]]) . Since it 
is a peaked function that rapidly declines to zero, it seems suitable 
to mimic what would be the expected behaviour of the metal de- 
position process. Nevertheless, it may be argued that the stochastic 
nature of the star formation algorithm, that samples the underlying 
"real" star formation on a discrete volume, should require weight- 
ing equally all the volume surrounding each stars particle. In order 
to judge the sensitivity of the resulting enrichment pattern on the 
choice of the weighting kernel, we will discuss in the following the 
effects of using instead a tof>-hat filter. 

Once the functional form of W is defined, one should choose 
the physical quantity with respect to which the weight is assigned. 
A first possibility is to weight according to the mass nig of each gas 
particle, in such a way that more massive particles receive a rela- 
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Figure 4. The effect of the numerical parameters, which define the distribution of metals, on global quantities computed within the virial radius of the CU 
cluster: total mass fraction of baryons in stars (left panel), mass-weighted Iron abundance (central panel) and fraction of metal mass in the diffuse gas (right 
panel). In the central and in the light panels, filled squares refer to Iron while open circles refer to Oxygen. The labels indicating the different runs are as 
reported in Table |2] 



tively larger amount of metals. An alternative choice is to weight 
instead according to the SPH volume carried by each particle, i.e. 
using the quantity nig/pg, where pg is the gas density associated 
to a particle. With this choice one gives more weight to those par- 
ticles which sample a larger volume and, therefore, collect a larger 
amount of metals, under the assumption that they are released in 
a isotropic way by each star particle. Therefore, one expects that 
weighting according to the volume assigns more metals to gas par- 
ticles which are relatively more distant from the star forming re- 
gions and, therefore, at lower density. 

In the following, we will assume in our reference run that the 
spreading of metals is performed with the spline kernel, by weight- 
ing over 64 neighbors according to the volume of each particle. 
In order to check the stability of the results, we also modified the 
scheme for metal and energy spreading in the following ways: (i) 
use the same kernel and density weighting, but with 16 and 128 
neighbours (N16 and N128 runs, respectively); (ii) weight accord- 
ing to the mass of the gas particle, instead of according to its vol- 
ume, using Nnb = 64 (MW run); (iii) use a top-hat window en- 
compassing 64 neighbors, weighting by mass (TH run). 

Figure |4] shows the global properties of the simulated cluster 
at z = 0, in terms of amount of stars produced, mass-weighted 
ICM metallicity and fraction of metals in the gas, for the different 
numerical schemes used to distribute metals. In general, we note 
that changing the details of the metal spreading has only a modest 
effect on the overall stellar population and level of enrichment. We 
find that the amount of stars within the cluster virial region ranges 
between 23 and 26 per cent of the total baryon budget. As for the 
global ICM enrichment, it is Zpe ~ (0.20 - O.25)ZFe,0 and 
Zo — 0.15^0,0 for the mass-weighted Iron and Oxygen abun- 
dance, respectively. Quite interestingly, only about a quarter of the 
total produced Iron is in the diffuse hot gas, while this fraction de- 
creases to about 15 per cent for Oxygen. This different behaviour 
of Oxygen and Iron can be explained on the ground of the differ- 
ent life-times of the stars which provide the dominant contribution 
to the production of these elements. Since Oxygen is produced by 
short-living stars, it is generally released in star forming regions 
and therefore likely to be distributed among star-forming gas parti- 
cles. For this reason, Oxygen has a larger probability to be promptly 
locked back into stars. On the other hand. Iron is largely contributed 
by long-living stars. Therefore, by the time it is released, the con- 



dition of star formation around the parent star particle may be sig- 
nificantly changed. This may happen both as a consequence of the 
suppression of star formation inside galaxies or because the star 
partic le may have later tak en part to the diffuse stellar component 
(e.g., iMurante et all2004l) . In both cases. Iron is more likely to be 
distributed among non star-forming particles and, therefore, con- 
tributes to the ICM enrichment, instead of being locked back in 
stars. 

This effect is clearly visible in the right panels of Figure [5] 
where we show the map of the fractional contribution of SNII to 
the global metal enrichment. Clearly, SNII provide a major contri- 
bution (magenta-red color) in high-density patchy regions, within 
and around the star forming sites. On the contrary, the contribution 
from SNIa and AGB dominates in the diffuse medium. This is a 
typical example of how the pattern of the ICM enrichment is de- 
termined by the competing effects of the chemical evolution model 
and of the complex dynamics taking place in the dense environment 
of galaxy clusters. It confirms that a detailed study of the chemical 
enrichment of the diffuse gas indeed requires a correct accounting 
of such dynamical effects. We find that the contribution in Oxygen 
within the cluster virial radius from SNII, from SNIa and from low- 
and intermediate mass stars is of about 70, 5 and 25 per cent respec- 
tively, while that in Iron is 25, 70 and 5 per cent. This demonstrates 
that none of these three sources of metals can be neglected in a de- 
tailed modelling of the chemical enrichment of the ICM. As shown 
in the left panels of Fig.|5] the distribution of Iron generally follows 
the global large-scale structure of gas inside and around the cluster, 
with an excess of enrichment inside the virial region and along the 
filaments from which the cluster accrete pre-enriched gas. We will 
comment in Sect. |3.3.T| the dependence of this enrichment pattern 
on the IMF. 

Quite interestingly, -Yam et al. I ( l2003l) discovered two 
SNIa not associated to cluster galaxies in Virgo and argued that 
up to about 30 per ce nt of the SNIa par ent stellar population is in- 
tergalactic (see also iMaoz et al]|2005l) . This is exactly the SNIa 
population that in our simulations is responsible for the diffuse en- 
richment in Iron. As the statistics of the observed population of in- 
tergalactic SNIa population improves, it will be interesting to com- 
pare it with the predictions of our simulations and to better quantify 
their contribution to the ICM enrichment. 

Figure |6] shows how the history of enrichment and of star for- 
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Figure 5. Maps of gas metallicity for t he three differe nt IMFs. Top, centeal and bottom panels refer to the ISalpeteJ jl955l) IMF, to the top-heavy 
iArimoto & Yoshi3ll987l) IMF and to the lKroup3 )200lh IMF, respectively. The left panels are for the Iron abundance, while the right panels show the 
fractional contribution of SNII to the global metallicity. For reference, the virial region of the cluster is marked with the white square. The maps have been 
done by smoothing the contribution of each gas particle with a SPH kernel using the corresponding smoothing length. The projection has been done through a 
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Figure 6. Left panel: the effect of changing the number of neighbours for metal spreading on the star formation rate (bottom) and on the SN rates (top; heavy 
and solid lines are for SNIa and SNII, respectively). Right panel: the enrichment history of Iron. The bottom panel shows the evolution of the mean Iron 
abundance in the gas; the top panel is for the evolution of the fraction of the total mass of Iron which is in the diffuse gas, the rest being locked in stars. 
In all panels, solid, dashed and dot-dashed curves are for the reference run (R), for the runs using 128 (N128) and 16 (N16) neighbors for metal spreading, 
respectively, when using the B-spline kernel with density weighting. 



mation changes by changing the scheme for distributing metals. 
In the left panel, the star-formation rate (SFR) is compared to the 
rate of SN explosions. As expected, the SNII rate follows quite 
closely the SFR (see left panel), and peaks at 2 ~ 2-3. On the 
other hand, the rate of SNIa is substantially delayed, since they 
arise from stars with longer life-times. Therefore, this rate peaks at 
« ~ 1.5 with a slower decline at low redshift than for SNII. Since 
the rate of SNIa is given by the combined action of SFR, life-time 
function and IMF, computing it by adding a constant time-delay 
to the SFR is generally too crude and inaccurate an approximation. 
As shown in Figurel?] the different redshift dependence of the SNIa 
and SNII rates is reflected by the different histories of enrichment 
in Iron and Oxygen. Since Oxygen is only contributed by short- 
living stars, the corresponding gas enrichment stops at 1, when 
the SFR also declines. On the contrary, the gas enrichment in Iron 
keeps increasing until 2 = 0, as a consequence of the significant 
SNIa rate at low redshift. An interesting implication of Fig|7]is that 
the relative abundance of Oxygen and Iron in the ICM can be dif- 
ferent from that implied by the stellar yields. Therefore, the com- 
monly adopted procedure to compare X-ra y measurements of rela- 
tive ICM abundances to stellar yields (e.g., iFinoguenov & PonmanI 
ll999l ; lBaumgartner et alj2"005h may lead to biased estimates of the 
relative contributions from SNIa and SNII (e.g., ?). We finally note 
that observations of SNIa rates in cluster s indicates fairly low rates, 
for both nearby and distant objects (e.g., iGal- Yam et al.l20"o3) . also 
consistent with the rates observed in the field. We postpone to a fu- 
ture analysis a detailed comparison between the observed and the 
simulated SNIa rates, and the implications on the parameters defin- 
ing the model of chemical enrichment. 

In general, using a different number of neighbours to spread 
metals has only a minor impact on the history of star formation and 
enrichment (see Fig. [6j. We only note that increasing the number 
of neighbors turns into a slight increase of the SFR at all redshifts, 
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Figure 7. The evolution of enrichment for the reference (R) run in both 
Iron (solid curves) and Oxygen (dashed curves). The bottom panel shows 
the evolution of the mean abundances in the gas, while the top panel is for 
the evolution of the fraction of the total metal mass which is in the diffuse 
gas, the rest being locked in stars. 



with a corresponding slight increase of the metallicity. In fact, in- 
creasing the number of neighbours has the effect of distributing 
metals more uniformly among gas particles. This causes a larger 
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Figure 8. The effect of changing the metal spreading on the Iron abundance profiles. Left panel: compaiison between the reference run (solid curve) and 
runs done by changing the number of neighbours to 16 (N16; short-dashed line) and to 128 (N128; long-dashed line). Right panel: comparison between the 
reference run (solid curve) and the runs in which metal spreading is performed using mass-weighting with the SPH kernel (MW; short-dash ed line) and the 
top-ha t kernel (TH; long-dashed line). The data points refer to a subset of four clusters, out of 16 observed with Chandra ana analysed by I Vikhlinin et alj 
12005 *). which have temperature in the range 2-A keV. For reasons of clarity, errorbars are not reported. 



number of particles to have a more efficient cooling and, therefore, 
to become eligible for star formation. 

In the left panel of Figure [S] we show the effect of changing 
the number of neighbors over which to distribute metals on the 
profile of the Iron abundance. As expected, increasing the number 
of neighbors corresponds to an increasing efficiency in distributing 
metals outside star forming regions. As a result, metallicity profiles 
becomes progressively shallower, with a decrease in the central re- 
gions and an increase in the outer regions. Although this effect is 
not large, it confirms the relevance of understanding the details of 
the mechanisms which determine the transport and diffusion of the 
metals. 

As for the comparison with data, we note that the differences 
between the different weighting schemes are generally smaller than 
the cluster-by-cluster variations of the observed abundance gradi- 
ents. In general, the simulated profiles are in reasonable agreement 
with the observed ones. A more detailed comparison with observed 
abundance profiles will be performed in a forthcoming paper, based 
on a larger set of simulated clusters (Fabjan et al. in preparation). 

Finally, we show in the right panel of Figure [8] the variation 
of the Iron abundance profile when changing the weighting scheme 
for the distribution of metals, while keeping fixed to 64 the number 
of neighbors. As for the Iron profile, using volume, instead of mass, 
in the SPH kernel has a rather small effect. Only in the innermost 
bin, the Iron abundance increases when weighting according to the 
mass as a result of the less effective spreading to less dense gas 
particles. As for the top-hat kernel, its Iron profile lies below the 
other ones at all radii, although by a rather small amount. 



3.2 The effect of resolution 

In this Section we present the results of the simulations of the C12 
cluster, done at three different resolutions (see Table[T). 

Figure[TO]shows the effect of resolution on the rates of star for- 
mation and SN explosions (left panel) and on the history of chemi- 
cal enrichment (right panel). As expected, increasing resolution en- 
hances the high-redshift tail of star formation, as a consequence of 
the larger number of resolved small halos which collapse first and 
within which gas cooling occurs efficiently. Quite interestingly, the 
increase with resolution of the high-z star formation rate is com- 
pensated by a corresponding decrease at lower, z"^ 1, redshift. As a 
net result, the total amount of stars formed within the cluster virial 
region by z = (left panel of Figure [9]( turns out to be almost 
independent of resolution. This results is in line with that already 
presented bv .Borgani et al. (2006) for a similar set of simulations, 
but not including the chemical enrichment. On the one hand, in- 
creasing resolution increase the cooling consumption of gas at high 
redshift, thereby leaving a smaller amount of gas for subsequent 
low-z star formation. On the other hand, smaller halos, forming at 
higher redshift when increasing resolution, generate winds which 
can more easily escape the shallow potential wells. As a result, the 
gas is pre-heated more efficiently, so as to reduce the later star for- 
mation. 

In spite of the stable star fraction, the level of ICM enrichment, 
both in Iron and in Oxygen (central panel of Fig. [9} increases with 
resolution. The reason for this is the larger fraction of metals which 
are found in the diffuse gas at increasing resolution (left panel of 
FiglSj. The fact that increasing resolution corresponds to a more 
efficient distribution of metals is also confirmed by the behaviour 
of the Iron abundance profiles (left panel of Figure fTTTl. which be- 
come systematically shallower at large radii, 0.5i?5oo • The rea- 
son for this more efficient spreading of metals from star-forming 
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The same as in Figure|4] but for the runs of the C12 cluster at three different resolutions. The labels indicating the different runs are as reported in 
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Figure 10. Left panel: the effect of resolution on the star formation rate (bottom) and supernova rates (top; tick lines: SNIa; thin lines: SNll). Right panel: the 
effect of resolution on the effect on the enrichment history. The upper panel shows the evolution of the fraction of Iron in contained in gas, while the lower 
panel shows the evolution of the Iron abundance in the gas. Solid, dashed and dot-dashed lines are for the Rl, R2 and R3 runs, respectively. 



regions has two different origins. First of all, by increasing reso- 
lution one better resolves processes, like "turbulent" gas motions 
and ram~pressure stripping, which are effective in diffusing met- 
als away from the star forming regions, thus preventing them to 
be locked back in stars. Furthermore, the better-resolved and more 
wide-spread star formation at high redshift releases a larger frac- 
tion of metals in relatively shallower potential wells. Therefore, 
galactic winds are more efficient in distributing metals in the IGM. 

As for the history of enrichment (right panel of Figure llOb , 
we note that increasing resolution has the effect of progressively 
increasing the gas Iron abundance at all redshifts. While the overall 
effect is of about ~ 30 per cent at 2 = it is by a factor of 2 or 
more at 4. This illustrates that, while resolution has a modest, 
though sizable, effect at low redshift, it must be increased by a sub- 
stantial factor to follow the enrichment process of the high-redshift 
IGM. As for the evolution of the fraction of Iron in gas, we note 



that it decreases with resolution at high redshift, while increasing at 
low redshift. The high-z behaviour is consistent with the increasing 
star-formation efficiency with resolution, which locks back to stars 
a larger fraction of metals. On the other hand, at low redshift this 
trend is inverted (see also the right panel of Fig[9f. This transition 
takes place at about the same redshift, z ~ 2.5, at which the star 
formation rate of the R3 run drops below that of the Rl runs, thus 
confirming the link between star-formation efficiency and locking 
of metals in stars. 

An increased efficiency with resolution in distributing metals 
in the diffuse medium is also confirmed by the the Iron abundance 
profile (left panel of Figure [TTt. While we do not detect any ob- 
vious trend at small radii, 0.2i?5oo, there is a clear trend for 
profiles to be become shallower at larger radii as resolution in- 
creases. In order to better show this effect, we plot in Figure [12] 
the Iron abundance profile out to 2i?vir, by using linear scales for 
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Figure 11. The effect of resolution on the abundance profiles of the C12 cluster. Left panel: profile of the mass-weighted Iron abundance, out to R^qq. The 
observational data points in the left panel are the same as in Figure[8] Right panel: the profiles of the relative abundance of Silicon (top) and Oxygen (bottom) 
with respect to Iron, plotted out to 2R^\^. The meaning of the different Une types is the same as in the left panel of Figure[TT] The dotted vertical line indicates 
the value of iisoo ■ 



the cluster-centric distance. This allows us to emphasize the regime 
where the transition from the ICM to the high-de nsity Warm-Hot 
Interg alactic Medium (WHIM) takes place (e.g., ICen & Ostrikeil 
I2OO6I) . Quite interestingly, the effect of resolution becomes more 
apparent in the outskirts of clusters. Behind the virial radius, the 
abundance of Iron increases by 50 per cent from the low-resolution 
(LR) to the high-resolution (HR) runs. In these regions the effect 
of ram-pressure stripping is expected to be less important, owing 
to the lower pressure of the hot gas. This demonstrates that the in- 
creasing of the ICM metallicity with resolution is mainly driven by 
a more efficient ubiquitous high-redshift enrichment, rather than 
by ram-pressure stripping of enriched gas from merging galaxies. 

As for the profiles of the relative abundance (right panel of 
Fig. II It. they are rather flat out to i?vir, with a relatively higher 
abundance for Silicon. In the innermost regions, the abundance ra- 
tios increase, with a more pronounced trend for [0/Fe]. The reason 
for this increase is an excess of recent star formation taking place at 
the cluster centre. As a consequence, elements mainly produced by 
short living stars, such as Oxygen, are released in excess with those, 
like Iron, which are mostly contributed by long-living stars. This 
also explains why the central increase is less apparent for [Si/Fe], 
being Silicon contributed by SNIa more than Oxygen is. An excess 
of star formation in the central regions of galaxy clusters is a well 
known problem of simulations, like those discussed here, which 
include only stellar feedback. For instance. ISaro et ^ 1 I2OO6I) anal- 
ysed a set of simulations, similar to those presented here, to study 
the properties of the galaxy population. They concluded that the 
brightest cluster galaxies (BCGs) are always much bluer than ob- 
served, as a consequence of the low efficiency of SN feedback to 
regulate overcooling in the most massive galaxies. 

In the outer regions, Rviv, the two relative abundances 
tend to increase, again more apparently for Oxygen. The outskirts 
of galaxy clusters have been enriched at relatively higher redshift 
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Figure 12. The effect of resolution on the profiles or Iron abundance, out to 
2i?vir ■ Linear scales are used for the radius, to emphasize the behaviour of 
the profiles in the cluster outskirts. 



(see also discussion in Sect. 13.3.31 below), when the potential wells 
are shallower and the enrichment pattern tends to be more uni- 
formly distributed. This causes the products of short-living stars to 
be more effectively distributed to the diffuse gas than at lower red- 
shift. In this sense, the increasing trend of the relative abundances 
produced by short- and long-living stars represents the imprint of 
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the different enrichment epochs. As for the dependence on resolu- 
tion, we note a systematic trend for an increase of Oxygen and, to 
a lesser extent, of Silicon, at least for 0. 17?5oo • This behaviour 
is consistent with the increased star formation rate at high redshift. 
Since Oxygen is relatively more contributed by short-living stars, 
then it is released at higher redshift than Iron and, therefore, has a 
more uniform distribution, an effect that increases with resolution. 
In the innermost cluster regions, _Ri$ O.I-R500, resolution acts in 
the direction of reducing the excess of Oxygen and Silicon. This 
can also be explained in terms of the dependence of the SFR on 
resolution: since most of the low-redshift SFR is concentrated at 
the cluster centre, its reduction at low redshift also reduces in these 
regions the relative amount of metals released from short-living 
stars. 

In conclusion, our resolution study demonstrates that the gen- 
eral pattern of the ICM chemical enrichment is rather stable in the 
central regions of clusters, O.Siisoo. However, the situation is 
quite different in the cluster outskirts, where both absolute and rel- 
ative abundances significantly change with resolution, as a conse- 
quence of the different efficiency with which high-redshift star for- 
mation is described. On the one hand, this lack of numerical con- 
vergence becomes quite apparent on scales -R500, which can be 
hardly probed by the present generation of X-ray telescopes. On 
the other hand, resolution clearly becomes important in the regime 
which is relevant for the study of the WHIM, which is one of the 
main scientific drive s for X-ray telescopes of the next generation 
(e.g., I Yoshikawa eral.,.2004,) . 



3.3 Changing the model of chemical evolution 

3. 3. 1 The effect of the IMF 

As already discussed in 12.3.41 a vivid debate exists in the liter- 
ature as to whether the level of the ICM enrichment can be ac- 
counted for by a standard, Salpeter-like, IMF or rather requires 
a top-heavier shape. The absolute level of enrichment in one el- 
ement, e.g. Iron, does not necessarily represent an unambiguous 
imprint of the IMF in our simulations. For instance, an exceedingly 
top-light IMF could still produce an acceptable Iron abundance in 
the presence of an excess of star formation in simulations. For this 
reason, it is generally believed that a more reliable signature of the 
IMF is provided by the relative abundances of elements which are 
produced by SNIa and SNII. 

As shown in Figure [Til the effect of assuming an IMF, which 
is top-heavier than the Salpeter one, is that of significantly increas- 
ing the number of SNII and, to a lesser extent, also the number of 
SNIa. This is consistent with the plot of Fig. [3] which shows that 
an Arimoto-Yoshii IMF predicts more stars than the Salpeter one 
already for M, ~ 1.5Mq. The larger number of SN clearly gener- 
ates a higher level of enrichment at all redshifts (bottom right panel 
of Fig. [14] A higher level of gas enrichment increases the cooling 
efficiency and, therefore, the star-formation rate (bottom left panel 
oflI4b. A higher star formation efficiency has, in turn, the effect of 
increasing the fraction of Iron which is locked in the stellar phase 
(top-right panel ofll4t. 

In Figure|5]we show the maps of Iron abundance (left panels) 
and of the fractional enrichment from SNII for the three IMFs. The 
effect of a top-heavy IMF is confirmed to increase the overall level 
of enrichment in Iron. At the same time, the contribution of SNII 
becomes more important, thus consistent with the increase of the 
number of massive stars. 

The effect of assuming a top-heavier IMF is quite apparent on 



the profiles of Iron abundance (left panel of FigurellSb. The level of 
enrichment increases quite significantly, up to factor of two or more 
at O.S-Rsoo, bringing it to a level in excess of the observed one. 
As expected, the larger fraction of core-collapse SN also impacts 
on the relative abundances (right panel of Fig. I15t . especially for 
[O/Fe]. Since Oxygen is largely contributed by SNII, its relative 
abundance to Iron increases by about 60 per cent. 

This result goes in the direction of alleviating the tension be- 
tween the largely sub-solar [O/Fe] values found for the Salpeter 
IMF and the nearly solar val ue reported by observational analy- 
ses (e.g., iTamura et al1l2004h . However, an overestimate of Oxy- 
gen from the spectral fitting, used in the analysis of observational 
data, may arise as a consequence of a temperature-dependent pat- 
tern of enrichment. Rasia et al. (in preparation) analysed mock 
XMM-Newton observations of simulated clusters, including chem- 
ical enrichment with the purpose of quantifying possible biases in 
the measurement of ICM metallicity. As for the Iron abundance, 
they found that its emission-weighted definition is a good proxy 
of the spectroscopic value. On the contrary, the spectroscopic mea- 
surement of the Oxygen abundance turns out to be significantly 
biased high with the respect to the intrinsic abundance. The reason 
for this bias is that, unlike Iron, the Oxygen abundance is obtained 
from emission lines which are in the soft part of the spectrum. On 
the other hand, relatively colder structures such as filaments, seen 
in projection and surrounding the ICM, give a significant contribu- 
tion to the soft tail of the spectrum. Since these structures are on 
average more enriched than the hot ICM, they are over-weighted 
when estimating element abundances from soft (;$ 1 keV) transi- 
tions. This is the case of Oxygen, whose abundance is generally 
estimated from the O-VIII line, which is at about 0.65 keV. 

Addressing the issue of observational biases in the measure- 
ment of the ICM enrichment is outside the scope of this paper. Still, 
the above example illustrates that such biases need definitely to be 
understood in detail if the enrichment pattern of the ICM has to be 
used as a fossil record of the past history of star formation in cluster 
galaxies. 

As for the results based on the IMF by iKroupal j200lh . we 
note that it induces only a small difference in the SFR with respect 
to the Salpeter IMF. While the SNIa rate is also left essentially un- 
changed, the SNII rate is now decreased by about 50 per cent. The 
reason for this is that the Kroupa IMF falls below the Salpeter one 
in the high mass end, ^ 5i\/0 (see Figure [3}. This is consistent 
with the maps shown in the bottom panels of Fig. |5] The global 
pattern of Iron distribution is quite similar to that provided by the 
Salpeter IMF, while the relative contribution from SNII is signifi- 
cantly reduced. Consistent with this picture, the profile of the Iron 
abundance shown in Fig. [TS] does not show an appreciable varia- 
tion with respect to the reference case. On the contrary, the profile 
of the Oxygen abundance and, to a lesser extent, of Silicon abun- 
dance, decreases significantly. 

Our results confirm that [a / Fe] relative abundances are sen- 
sitive probes of the IMF. However, we have also shown that dif- 
ferent elements are spread in the ICM with different efficiencies. 
This is due to the fact that long-lived stars can release metals away 
from star forming regions and, therefore, their products have an en- 
hanced probability to remain in the diffuse medium (see discussion 
in Sect. 13. II ). Therefore, both a correct numerical modeling of such 
processes and an understanding of observational biases are required 
for a correct interpretation of observed metal content of the ICM. 
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Figure 13. The same as in Figure|4] but changing parameters related to the model of chemical evolution and to the feedback. The labels indicating the different 
runs are as reported in Table |2] 
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Figure 14. The same as Figure[6] but cha nging t he IMF. The sol id lines correspond to the reference run, while the short-dashed and the long-dashed lines are 
for the IMFs bv lArimoto & Yoshiil h987l) and bv lKroupj feOOlh . respectively. 



3.3.2 Changing the yields 



In Figure 1161 we show the effect of using the yields by 
Chiefs & Limongil l l2004h (CL) for the SNII, instead of those by 
Wooslev & Weaved ( Il99a) (WW) as in the reference run, on the 
Iron density profile. As for the enrichment pattern of the ICM, using 
either one of the two sets of yields gives quite similar results, both 
for the abundance of Iron and for the [0/Fe] relative abundance. 
This is apparently at variance with the results shown in FigH] where 
we have shown the differences in the production of different metals 
from an SSP for the two sets of yields. However, we note from that 
figure that such differences in the metal production have a non triv- 
ial dependence on the SSP characteristics, being different elements 
over- or under-produced for one set of yields, depending on the 
SSP initial metallicity. On the other hand, the metallicity of each 
star particle in the simulations depends on the redshift at which it 



is created, with the star-formation history being in turn affected by 
the enrichment pattern. 

As for the enrichment in Iron (left panel of Fig.|16t, we remind 
that this element is mostly contributed by SNIa, while the contri- 
bution from SNII is not only sub-dominant, but also preferentially 
locked back in stars. Since we are here changing the yields of the 
SNII, there is not much surprise that the effect of the profiles of 

is marginal. The situation is in principle different for Oxygen 
(right panel of Fig|16t. However, Fig|2] shows that the WW yields 
for Oxygen are in excess or in defect with respect to the CL ones, 
depending on the initial metallicity. As a result, we find that [0/Fe] 
for the diffuse gas is left again substantially unchanged, while a sig- 
nificant variation is found for stars, whose [0/Fe] is about 40 per 
cent larger when using the CL yields table. The profiles of [0/Fe] 
for stars show a mild decrease at large radii, consistent with the fact 
that Oxygen is more efficiently spread around stars in the regions 
where enrichment takes place at higher redshift (see discussion in 
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Figure 15. The effect of changing the IMF and the yields on the abundance profiles. Left panel: profiles of the mass-weighted Iron abundance. The data points 
are the same as in Figure[8] Right panel: the relative abundance of Silicon (top) and Oxygen (bottom) with respect to Iron, for the same runs shown in the left 
panel. The different line-types have the same meaning as in Figure [T4l 



Sect. 13. lb . Once again, this difference in the enrichment pattern of 
gas and stars reflects the different efficiency that gas dynamical pro- 
cesses have in transporting metals away from star forming regions. 

As a final remark, we would like to emphasize that several 
other tables of stellar yields have been presented in the literature, 
besides the ones by WW and C L considered here, which refer 
only t o massiv e stars, and those b y van den Hoek & GroenewegenI 



only t o massiv e stars, and those by Ivan den hloek & uroenewegenl 
( Il997 l) and by iThielemann et al.l HoOSh that we adopted for low 
and intermediate mass stars and for SNla, respec tively. Besides 
the y ields for intermediate-mass stars computed by Renzini & Volil 
( 1198 li) , other sets of metallici ty-d ependen t yield s have been pro- 
vided bv llwamoto et alj i 19991) andlMarigo j200lh for intermediate 



and low-mass stars, and bv'Po rtinari et al.l d 1998 1) for massive stars. 
The differences among such sets of yields are due to the different 
stellar evolutionary tracks used and/or to the different way of de- 
scribing the structure of the progenitor star . A different approach 
has been followed by [Francois et alj ( |2004) . who inferred stellar 
yields by requiring that their model of chemical evolution repro- 
duces the observed enrichment pattern of the Milky Way. 

Assessing in detail which one of them should be preferred to 
trace the cosmic evolution of metal production with hydrodynam- 
ical simulations is beyond the scopes of this paper. It is however 
clear that any uncertainty in the models on which yields computa- 
tions are based can pr oduce significant c hanges in the pattern the 
ICM enrichment (e.g.. lGibson et al.lll997h . 



3.3.3 The effect of the lifetimes 

As we have sh own in Figure [T] using the lifetime function by 
iMaeder & Me vnel ( 1989), instead of that by Padova ni & Matteu cci 
1 19931) , corresponds to significantly increasing the delay time 
for the release of energy and metals from low-mass stars with 
M'^ IMq, while leaving substantially unchaged that of massive 
stars. As a result, we expect that the SNII rate will remain substan- 



tially unchanged, while shifting towards lower redshift the rate of 
SNIa. This expectation is completely in line with the results shown 
in the left panel of Figure (TT] changing the lifetime function has a 
small enough effect on the history of enrichment that the resulting 
SFR and SNII rate are left completely unchanged; on the other hand 
the SNIa rate is suppressed at high redshift, while it is substantially 
increased at 0.5. 

The consequence of changing the lifetimes on the typical 
epoch of the ICM enrichment can be can be appreciated from Fig. 
[Tsl To this purpose, we computed for each gas particle the mean 
value of the cosmic time at which it has been enriched. This quan- 
tity is obtained by weighting each epoch, at which an enrichment 
episode takes place, by the mass of metals received within that 
time-step. In Figure [TS] we plot the radial profile of the mean cos- 
mic epoch of enrichment, which is computed by averaging over 
all the gas particles falling inside different radial intervals. Outside 
the virial radius, for both lifetime functions the enrichment is pro- 
gressively more pristine as we go to larger cluster-centric distances. 
While at 7?virthe mean redshift of enrichment is about 0.7, this in- 
creases to ^ 1.5 at 2i?vir. Quite interestingly, inside ~ 0.8J?virthe 
age of enrichment becomes nearly constant down to the central re- 
gions, ~ 0.2i?vir. This change of behaviour further indicates that 
processes other than star formation, e.g. related to gas dynamical 
processes, become important in determining the enrichment pat- 
tern. In the core region, ^ 0.2i?vir, there is a sudden decrease of 
the enrichment redshift. This is due to the excess of recent star for- 
mation, which takes place in the central cluster galaxy, and which 
also causes the spike in the [0/Fe] relative abundance. 

Using the lifetimes bv lMaeder & Mevnetl ( [l98^ turns into a 
more recent enrichment of the ICM. We should remember here that 
this figure shown the mean age at which gas has been enriched in 
all metals. Since the global metallicity is dominated by Oxygen, 
mostly produced by SNII, we expect a more pronunced effect on 
the enrichment in Iron. In principle, a more recent enrichment due 
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Figure 16. The effect of changing the SNII yields on the abundance profiles. Left panel: profile of the Fe abundance in the ICM; the data points are the same 
as in Figure[8] Right panel: the profiles of the [0/Fe ] relative abundance s for th e ICM (lower panel) and for the stars (up per panel). Solid aiid dash ed curves 
refer to the reference (R) run, based on the yields bv lWooslev & WeaveJ fT99^ . and to a run (CL) in which the yields bv lChiefii & Limond j2004l) are used 
instead. 




Figure 17. Left panel: the effect of changing the lifetime function on the star formation rate (bottom panel) and supernova rates (top panel; tick lines: SNIa; 
thin lines: SNII). The meaning of the different line types is the same as in the left panel of Figure [TT] Right panel: the effect of the lifetime function on the 
profile of the Iron abundance (right panel). The data points are the same as in Figure [8] The solid fine is for the reference (R) run, which uses the lifetime 
function bv .Padovani & Matteucci tl993,) . while the dashed line is for the MM run which uses the lifetime function bv .Maeder & Mevnet tl989.) . 



to longer delay times of SNIa could help explaining the observa- 
tional evid ences for an evolution of the ICM iron content at 0.5 
tealestra et al. 2007). We will postpone a more detailed compari- 
son with the observed ICM metallicity evolution to a forthcoming 
paper (Fabjan et al., in preparation). 



As for the profile of the Iron abundance, we note that there 
is a small but sizeable increase, especially at small radii. This in- 
crease is due to a more efficient unlocking of Iron from star form- 
ing regions. Since the delay time is enhanced for SNIa, they have 
an larger probability to release metals outside star forming regions. 
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Figure 18. The radial dependence of the mean value of the age of the Uni- 
verse at which the ICM has been enriched. The solid curve is for the ref - 
erence (R) runs, which uses the lifetimes bv lPadovani & MatteucciHl993h. 
while the dashed curve is when using the lifetimes bv iMaeder&Mevnell 
Horizontal dashed lines indicate different values of redshift corre- 
sponding to difference cosmic times. 

On the other hand, we verified that the enrichment is Oxygen is left 
unchanged, as expected from the stability of the SNII rates. As a 
result, we find that [O/Fe] decreases by about 0. 1 at i?5oo and by 
0.2 in the innermost regions. 



efficiency also provides a more efficient distribution of metals in 
the gas. This is confirmed by the maps of Iron abundance shown in 
Figure [21] In the absence of winds, the pattern of the Iron distri- 
bution is quite patchy, with lumps of high abundance surrounding 
the star-forming regions. On the other hand, including the effect 
of winds turns into a more diffused enrichment pattern, an effect 
which increases with the winds' speed. 

As shown in the right panel of Fig. [191 this effect is more pro- 
nounced for Oxygen, whose mass fraction in gas increases by a fac- 
tor of about 2.5, than for Iron, for which the same fraction increases 
by about 50 per cent. As already discussed. Iron is more efficiently 
distributed outside star-forming regions. This is due to the fact that 
Iron is largely contributed by long-lived stars, which explode when 
star particles possibly had the time to leave star-forming regions. 
Therefore, these regions are expected to be relatively more enriched 
in Oxygen. Since galactic winds are uploaded with gas from star- 
forming particles, they are relatively more enriched in Oxygen. As 
a consequence, increasing the feedback from winds has the general 
effect of increasing the diffusion of metals, by an amount which is 
larger for Oxygen than for Iron. This effect also explains why the 
ICM metallicity is not suppressed in the same proportion of star 
formation, when feedback from winds increases (see left and cen- 
tral panel of Fig.[T9ll. 

As shown in the right panel of Fig|20] the effect of winds 
in unlocking metals from stars is more effective at high redshift. 
This is quite expected, owing to the shallower potential wells 
which are present at earlier times and make winds more effec- 
tive in transporting metals outside star-forming regions. This re- 
sult further illustrate the fundamental role that winds have played 
in establishing the enrichment pattern of the high-redshift IGM 
(iDave & Oppenheimer.,2006) . 



3.3.4 Changing the feedback strength 

Feedback from galactic winds has been originally implemented 
in the GADGET code with th e purpose of preventing overcooling 
dSpringel & Hernquislll20033) in simulations of galaxy formation. 
Besides producing a realistic star form ation in cosmological simu- 
lations jSpringel & Hernauislll2003bl) . galactic winds are also ex- 
pected to play a crucial role in the transport of metals outside 
star-forming regions and, therefore, in the enrichment of the dif- 
fuse gas. According to the implementation of galactic ejecta in 
the GADGET code, gas particles to be uploaded in winds are cho- 
sen among star forming particles which, by definition, are those 
more metal rich. Recent numerical studies of the key role played 
by winds in enriching diffuse co smic baryons have been presented 
bv lOave & Oppenheimej j2006t), who focused their analysis on the 
high-redshift IGM, and bv ICen & Ostrikeil ( |2006|) . who concen- 
trated on the study of the WHIM. 

The left panel of Figure [19] clearly shows that neglecting the 
effect of winds increases the star fraction in clusters up to ~ 35 
per cent. This value is possibly even an underestimate, owing to the 
lack o f numerical converge nce of the star fraction in the absence of 
winds lBorgani et al] ( 12006!) . On the other hand, increasing the wind 
speed to 1000 kms~^ reduces the fraction of cooled gas by almost 
a factor two. This reduction of the star fraction inside the cluster 
virial region at z = is just the consequence of the suppressed 
star formation due to the wind action (see the left panel of Figure 
I20t . As shown in the right panel of Figure [19] an increasing wind 



The dependence of the abundance profiles on the feedback 
strength is shown in Figure [22] The effect of increasing the wind 
speed from 500kms~^ to 1000 kms^^ is that of suppressing the 
profile of Iron abundance, as a consequence of the suppressed star 
formation, bringing it below the level indicated by observational 
data. On the other hand, the effect of neglecting galactic winds has 
a less direct interpretation. While for Rf^ 0.3Ji5oo the enrichment 
level is higher, due to the larger amount of stars, it falls below the 
profile of the reference (R) run at smaller radii, due to the reduced 
efficiency of metal transport. Only in the innermost regions, the 
excess of star formation, taking place in the absence of efficient 
feedback, causes a spike in the Iron abundance. 

As for the profiles of the relative abundances of [O/Fe] and 
[Si/Fe] (right panel of Fig|22t. they confirm the different effect that 
feedback has in determining the enrichment level for different el- 
ements. Consistent with what shown in Fig[T9] Silicon and Oxy- 
gen becomes progressively over-abundant with respect to Iron as 
the feedback strength increases. In the innermost regions, the ex- 
cess of low-redshift star formation in the absence of efficient feed- 
back manifests itself through a sudden increase of the relative abun- 
dances. A comparison of this figure with the right panel of Fig|15l 
shows that the effect of a stronger feedback is similar to that in- 
duced by assuming a top-heavier IMF. This highlights once more 
that using the observational determination of the [a/Fe] abundances 
to infer the IMF requires a good knowledge of other gas-dynamical 
effects. 
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Figure 19. The same as in Figure|4] but changing the velocity of the galactic winds. The labels indicating the different runs are as reported in Table[2] 
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Figure 20. The same as Figure|6] but changing the feedback strength associated to galactic winds. The solid lines correspond to the reference run, while the 
short-dashed and the long-dashed lines are for the runs with strong winds (SW) and switching off winds (NW), respectively. 



4 DISCUSSION AND CONCLUSIONS 

We presented results from hydrodynamical simulations of galaxy 
clusters, aimed at studying the metal enrichment of the Intra- 
Cluster Medium (ICM). Our simulations were based on a version 
of th e Tree-SPH GADGET-2 code ( Springel et al. 2001; Springel 
I2OOS ). which includes a detailed treatment of the production of 
metals from different stellar populations. In particular, we account 
for the delay times at which different stellar populations release 
energy and metals from SNIa and SNII explosions and stellar mass 
loss. The resulting chemo-dynamical code is rather flexible since 
it allows one to change the assumed stellar initial mass function 
(IMF), the life-time function and the tables of stellar yields. Thanks 
to the efficient implementation of the way in which the equations of 
the chemical evolution model are solved, the corresponding com- 
putational cost is only of about 10 per cent of the total cost of a 
simulation in a typical configuration. 

The analyses presented in this paper have been carried out 



with the twofold aim of assessing in detail the impact of numer- 
ical effects and of changing the parameters defining the chem- 
ical evolution and the feedback model on the resulting pattern 
of ICM enrichment. For this reason, we have concentrated our 
analysis only on two medium poor clusters, having virial masses 
Afvir — 2 X lO^** h~^MQ, which have been simulated in a vari- 
ety of conditions. For this reason, we have preferred not to perform 
here any detailed comparison with observational X-ray data. Such 
a comparison would require simulations of a more extended set 
of clusters, so as to appreciate both the cluster-to-cluster variation 
of the enrichment pattern and the dependence on the cluster mass 
(Fabjan et al., in preparation). 

The main results of our analysis can be summarized as fol- 
lows. 

(a) Gas-dynamical effects, related to gas mixing, galactic winds 
and ram-pressure stripping, play an important role in determin- 
ing the distribution and relative abundances of different chemical 
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Figure 21. Maps of the Iron metallicity for the run with no winds (NW, left panel), for the reference run (R, central panel; = 500 km s ^ ) and for the run 
with strong winds (SW, right panel; Vw = 1000 km s""^). Details of the maps are the same as in Fig|5] 
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Figure 22. The effect of changing the feedback strength on the abundance profiles. The dot-dashed lines are for the run with strong winds (SW), while the 
dashes lines are for the run without winds (NW). Left panel: profiles of the mass-weighted Iron abundance. The data points are the same as in Figure[8] Right 
panel: the relative abundance of Silicon (top) and Oxygen (bottom) with respect to Iron, for the same runs shown in the left panel. 



species. For instance, we find that the products of SNII are gener- 
ally more concentrated around star-forming regions than those of 
low-mass stars, due to the different time-scales over which stars 
of different mass release metals. This emphasizes the relevance of 
properly including in the simulations a correct description of the 
lifetimes of different stellar populations. 

(b) A compar i son w ith observed metallicity profiles obtained by 
IVikhlinin et"^ feOOSh from Chandra data shows a reasonable and 
encouraging agreement in both shape and normalization, when us- 
ing a Salpeter ( 1955) IMF, within ~ 0.5i?vir. The [0/Fe] relative 
abundance is predicted to be undersolar, ~ —0.35 over the cluster 
virial region, with a mild increase at larger radii and with a sudden 
increase in the core region, due to an excess of recent star forma- 
tion. 



(c) For each gas particle we estimated the metallicity-weighted 
cosmic time of enrichment. We find that this typical age of enrich- 
ment of the hot gas corresponds to redshift z ~ 0.5, over the whole 
cluster virial region. Only in the innermost region, which is heavily 
affected by the presence of the Brightest Cluster Galaxy (BCG), the 
enrichment took place at a significantly more recent epoch, due to 
an excess of low-redshift star formation. 

(d) As expected, the IMF has a strong impact on the enrichment 
level of both the ICM and the stellar populations. Assuming a tog- 
heavy IMF ( Ai'imoto & Yoshii 1987), instead of the Salpeter ( 1953) 
one, it turns into an increase of the Iron abundance and of the [0/Fe] 
relative abundance. Within the cluster virial regions, the Iron metal- 
licity of the ICM increases by almost a factor of two, while Oxy- 
gen increases by about a factor of three. Using the IMF proposed 
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bv lKroupj l200lh . it turns into a general decrease of [0/Fe] with 
respect to the prediction of a Salpeter IMF, owing to the relatively 
smaller number of massive stars, 
(e) Increasing the velocity of galactic ejecta from 500 kms^^ to 
1000 km s^^ provides a strong suppression of the star formation 
history and of the corresponding level of enrichment. The fraction 
of baryonic mass in stars within i?vir drops by about a factor of 
two, while the level of Iron enrichment drops only by about 40 per 
cent. This is due to the increasing efficiency of stronger winds to 
unlock metals from star-forming regions. 

In order to judge the stability of our results against numer- 
ical effects, we have considered several different prescriptions to 
distribute metals around stars and the effect of progressively in- 
creasing resolution. For instance, we find that in the central regions, 

0.3i?5oo, the pattern of ICM enrichment is quite insensitive to 
the shape of the weighting kernel, to the number of neighbours over 
which metals are distributed. The only significant trend found is for 
the profiles of the Iron abundance to become higher at larger radii 
as the number of neighbors used for the spreading increases. As a 
general result, we can then conclude that our results are not strongly 
affected by the details with which metals are distributed around 
star-forming regions. Furthermore, we find no obvious trend in the 
Iron abundance profile with resolution within 0.3i?5oo- However, 
we do find a significant and systematic increase of the enrichment 
level with resolution in the cluster outskirts: at Rft i?vir the Iron 
abundance in the highest resolution run is about 60 per cent higher 
that in the low resolution run. These differences are not due to a 
lack of convergence of resolution of the total amount of stars found 
&t z = Q. Instead, they are mainly due to the better-resolved star 
formation at high redshift, which provides a more uniform enrich- 
ment of the inter-galactic medium (IGM). This calls for the need 
of significantly increasing the resolution to have numerically con- 
verged results on the enrichment pattern in the cluster outskirts and 
in the Warm-Hot Intergalactic Medium (WHIM). 

The results presented in this paper and those from previous 
chemo-dynamical models of the ICM (e.g., see a lso iLia et al.l 
I2OO2I : IValdaminill2003l : [s ommer- Larsenetalj|2005h demonstrate 
that sophisticated models of chemical evolution can be efficiently 
included in advanced codes for cosmological simulations of galaxy 
clusters. This opens the possibility of accurately describing the his- 
tory of metal production through cosmic epochs, in a variety of 
environments. While we have focused our analysis on the low- 
redshift enrichment of the ICM, a number of other studies can 
be foreseen, from the enrichment of the IGM at high redshift, 
2 , to the metal abundances in both elliptical and disk galaxies 
(e.g. iKawata & Gibsoiill2003l : lKobavashi|[2004 IScannapieco et al.l 
l2006h . 

Even restricting to the study of the ICM, a number of direc- 
tions of further developments and improvements with respect to the 
analysis presented here can be devised. 

As for the comparison with observational data, a fundamen- 
tal step is represented by understanding in detail possible obser- 
vational biases in the measurement of the ICM metallicity. In the 
analysis of observational data, the ICM metallicity is estimated 
by fitting the X-ray spectrum to a plasma model. On the other 
hand, in the analysis of simulations, metallicity is usually esti- 
mated by either mass-weighting or emission-weighting the metal- 
licity carried by each gas particle. The best way of approach- 
ing this problem is by performing mock observations of simu- 
lated clusters, which reproduce as close as possible the observa- 
tional setup (i.e., instrument response and PSF, instrumental back- 



ground, etc.; Rasia et al., in preparation), and analyse them exactly 
in the same way as observational data. The preliminary compari- 
son with observational data, presented in this paper, shows an en- 
couraging agreement with the profiles of the Iron abundance as ob- 
taine d from Chandra obse rvations of a set of nearby relaxed clus- 
ters dVikhlinin et al.ll2005h . As a word of caution, we point out that 
for a numerical model to be fully successful, it should reproduce 
at the same time the properties of the ICM and those of the cluster 
galaxy population. However, a well known problem of simulations 
based on stellar feedback, like those presented here, is that they 
produce central cluster galaxies which are much bluer and more 
star-forming than observed (e.g., ISaro et al . 2006). In this respect, 
AGN feedback is generally considered as the most credible candi- 
date to quench star formation at low redshift, thereby accounting 
for the observed properties of both the brightest cluster galaxies 
and of the "cool cores". Besides regulating star formation, AGN 
feedback should also play an important role in distributing met- 
als from the core regions, through the generation of buoyant bub- 
bles created by collimated jets, which shock on the dense ambi- 
ent gas (e.g.. iDalla Vecchia et al.ll2004 ISiiacki & Springelll20063 : 
iRoediger et alj2007l . and references therein). 

From the observational side, the sensitivity of the present gen- 
eration of X-ray satellites are allowing us to measure the distri- 
bution of metals at most out to about i?5oo for nearby, z ~ 0.1, 
clusters, while global measurement of the IC M metallicity are pos - 
sible for distant cluster, out to 2 ~ 1 (e.g., ISalestra et al.ll20()7h . 
X-ray telescopes of the next generation, thanks to a much larger 
collecting area, improved spectral resolution and lower instrumen- 
tal background, will open a new window on our knowledge of the 
chemo-dynamical history of the cosmic baryons. Deep pointings 
of distant clusters will provide the signature of the metal enrich- 
ment in the just-assembling ICM, out to z ~ 2, thus bridging the 
gap between observations of the metal-enrichment of the low-z 
ICM and of the high-z IGM. There is little doubt that numerical 
simulations, like those presented here, will provide the ideal tool 
for the interpretation of these observations within the framework of 
cosmological models of structure formation. 
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